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ABSTRACT 


Basic  ideas  of  linear  system  theory  are  applied  to 
acoustic  scattering  problems.  A transfer  function 
and  an  impulse  response  for  general  scattering  bodies 
are  found  both  analytic  silly  and  experimentally.  The 
relationship  between  target  strength  and  the  transfer 
function  is  considered.  Explosive  pulses  are  dis- 


cussed briefly. 
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by 

David  L.  Clifton 


I.  LINEAR  SYSTEM  THEORY 


Definition: 

A linear  system  is  a set  of  specifications  for  determining  a 
function  y(t)  (output)  from  a given  function  x(t)  (input). 

y(t)  * L[x(t)]  , (1) 

where  L is  a linear  operator. 

Definition: 

If  L does  not  change  in  time,  and  if  the  input  is  a unit 
impulse,  then  the  output  is  known  as  the  impulse  response  of  the 
system. 

If 

L / L(t)  , and 

x(t)  = 6(t)  , 

then 

L[x(t)]  = y(t)  s h(t)  , (2) 

where  h(t)  is  the  impulse  response. 
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If  the  impulse  response  h(t)  of  a linear  system  is  known,  the  output 
y(t)  for  a given  input  x(t)  can  always  be  found.  The  relation  is 


x(T)h(t 


x(t  - t)h(T)dT 


(5 


which  is  often  written  as 

y(t)  = x(t)  * h(t) 

The  asterisk  is  read  "convolved  with". 

Equation  (3)  follows  from  the  property  of  the  delta  function 

that 


x(x)6(t  - r)dx 


Operating  on  both  sides  of  the  above  with  the  linear  operator  L,  and 
noting  that  L[x(t)]  = y(t)  and  L[6(a)]  = h(a),  the  result  is 


x(x)h(t  - x)dx 


Definition: 

Causal  system — In  a physical  passive  system,  the  output  is 
zero  for  t < a whenever  the  input  is  zero  for  t < a.  Such  a 
system  is  said  to  be  "causal". 


In  a causal  system,  the  impulse  response  h(t)  is  zero  for  t < 0. 
With  this  in  mind,  Eq.  (3)  ma y be  rewritten 


then 


PROOF 


Taking  the  Fourier  transfom  of  both  sides 


F[f(t)] 


F[f(t)] 


P[f(t)] 


’ 


Now 


e-la>(t-f) 


f2(t)dt  = Hf2(t)] 


Thus 


F[f(t)]  = f°°  e~^  f1(T)F[f2(t)]dT 

U-  00 


Definition: 

System  function  or  transfer  function- -that  function  which,  when 
multiplied  by  the  transform  (Fourier,  Laplace)  of  the  input,  yields 
the  transform  of  the  output. 

Given  a linear  system  with  input  x(t),  output  y(t),  and 
impulse  response  h(t). 

Fran  Eq.  (3) 


x(T)h(t  - x)dT 


Taking  the  Fourier  transform  of  both  sides, 


F[y(t)]  = F[x(t)]F[h(t)] 


(T) 


From  the  definition,  it  is  clear  that  F[h(t)]  is  the  system  function 
of  the  linear  system 

H(cu)  b system  function  = F[h(t)]  . (8) 
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Note  that  when  the  transfer  function  is  known,  the  output 
corresponding  to  a given  input  can  always  be  found.  The  relation 
between  the  output  of  a given  linear  system,  its  input,  and  its 
transfer  function,  is 

y(t)  = F _1^F{x(t)]HHj  . (9) 

Definition: 

A coherent  system  obeys  the  relation 

m 

F[y(t)]  = H(cu)F[x(t)]  , (lO) 

whereas  an  incoherent  system  obeys  the  relation 

|F[y(t)]|2  = |H(u>)|2  |F[x(t)]|2  . (11) 

Definition: 

The  average  power  of  a signal  x(t)  is  defined  by  the  limit 

*"2(t)  = ^ ^ l*(t)|2dt  . (12) 

Definition: 

For  a signal  with  zero  average  power,  hence  finite  energy, 
the  autocorrelation  function  is  defined  as 


x(t)x(t 


t 


00 

x(t)x(t  + t)dT  = R (-t) 

oo 


Definition: 

For  signals  with  zero  average  power,  the  cross- correlation 
functions  are 
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Henceforth  only  signals  with  zero  average  power  will  be  considered 


Consider  a linear  system  with  impulse  response  h(t) 
input  x(t),  and  output  y(t).  It  will  be  shown  that 


PROOF 


By  definition 


From  Eq.  (3) 


Substituting  this  into  the  above 


x( x - a)h(  c)dodT 


then 


x(p  + t - ff)h(  <x)dodP 


i < > 

Interchanging  the  order  of  integration, 


R (t)  = [ h(cr)  f x(p)x[p  + (t  - <r)dpda 
Y*  J- 00  J-oo 


From  the  definition  of  the  autocorrelation  function,  this  can  be 
rewritten  as 

/OO 

^(■t  - a)h(a)dc  , or 

CO 

Eyx(t)  ’ • 

Definition: 

The  power  spectral  density  Sxx(oo)  of  a function  x(t)  is 
defined  as  the  Fourier  transform  of  its  autocorrelation  function. 

sx>>  « . (16) 

Consequently,  the  autocorrelation  function  is  the  inverse  transform 
of  the  power  spectral  density. 

Definition: 

The  cross  power  spectral  densities  (CPSD)  of  two  functions 
x(t)  and  y(t)  are  defined  as  the  Fourier  transforms  of  their 
cross- correlation  functions . 


7 


j 

J 


a 


a 


rev  ..vfUL*- 


dt 


, and 


(17) 


S (os ) = f R (t)e_:kiyt  dt  . (18) 

xyv  ' J_M  xyv  ' K ' 


Thus  the  cross-correlation  functions  are  the  inverse  transforms  of 
the  CPSD's. 

Definition: 

White  noise  is  a signal  having  a power  spectral  density  with 
constant  amplitude  and  randomly  varying  phase . 

Consider  the  linear  system  characterized  by  input  x(t), 
output  y(t),  and  impulse  response  h(t).  It  will  be  shown  that  if 
x(t)  is  white  noise, 


where 


aR  (t)  = h(t) 
yx  w 

a = constant 


(19) 


PROOF: 

Rewriting  Eq.  (15), 


V(t>  - E**(t)  * h(t)  • 

Taking  the  Fourier  transform  of  both  sides  of  Eq.  (15),  making  use 
of  the  rule  established  in  Eq.  (6), 

SyxM  = Sxx(o>)H(o>)  . (20) 

If  x(t)  is  white  noise,  S is  a constant.  (Let  S (u>)  = — .) 

XX  XX  EL 
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For  such  an  input. 


“Syxfa)  = H(o>)  . (21) 

Taking  the  inverse  transform  of  both  sides  of  Eq.  (21), 

h(t)  = aRyx(t) 


II . ANALYTICAL  DETERMINATION  OF  THE  TRANSFER  FUNCTION 


l 


Consider  an  underwater  object  of  arbitrary  shape  with  a 
surface  S enclosing  the  origin.  Let  a pulse  x(r,t)  be  incident 
on  the  object  (along  the  line  BO),  and  designate  the  scattered 
pulse  by  y(r,t).  (See  Fig.  1.)  Assume  x(r,t)  is  composed  of  plane 
waves  moving  in  the  -a  direction,  and  designate  its  Fourier  trans- 
form by  X(r,co).  The  transform  of  y(r,t)  is  Y(r,o)).  Let  r^  indicate 
the  position  at  which  the  incident  and  scattered  pulses  are  observed. 
The  transfer  function  for  the  scattering  process,  observed  at  r^  is 


H(r  ,oj)  = 


Y(yq>) 

X(r  ,o>) 


(22) 


Given  the  incident  signal  x(r,t),  X(r^,cn)  can  be  easily 
obtained.  Y(r ,o>)  can  be  obtained  as  follows.  x(r,t)  is  broken 

"G 

into  its  component  sinusoids,  one  of  which  is 

P(r,  t)  = X(r,03)eKDt  ; (25) 

The  scattered  field  p(r,t)  due  to  this  steady- state  sinusoid  is 
found  by  solving  the  wave  equation. 


p(r,t)  = g(r,oj)elxut 

The  scattered  pulse  is  found  by  integrating  Eq.  (24)  over  all 
frequencies . 


(24) 


y(r,t)  =J 


\ loot  , 

g(r,cn)e  da> 


(25) 
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Without  carrying  out  the  integration,  it  will  be  possible  to  pick 
Y(r,u>)  out  of  the  integral  expression  for  y(r,t).  Thus  in  this  case, 

Y(r,oj)  = 2«  g(  r,u))  . (26) 

Knowing  Y(rt,cs)  and  X(r^,cu),  H(rt,o>)  can  be  found  from  Eq.  (22). 

Note:  The  position  of  the  observer  with  respect  to  the  target  (r  ) 

is  assumed  a constant  parameter  of  the  scattering  system.  Though 
it  is  indicated  as  a variable,  it  is  not  so  treated  in  the  calcu- 
lations . 


A.  Example  1:  Transfer  Function  for  a Soft  Cylinder 


Consider  an  infinite  soft  cylinder  of  radius  a,  with  its 
axis  along  the  y-axis.  (See  Fig.  2.)  The  incident  pulse  if  moving 
in  the  -Q : direction,  and  the  receiving  transducer  is  located  at  the 
point  Q.  Assume  an  incident  pulse  of  the  form 


ioi  t its  — 

x(r,t)  = A(ct  + orx)  e 0 e ° C , (27) 

where 

A(ct  + oP" ) = 1 if  |ct  + | S j 

A(ct  + a1)  = 0 if  |ct  + a1 1 > l , 

1 

a = or  + b , and 


b = constant 


The  Fourier  transform  of  x(r,t)  is 


. — sin(o3  - us  ) 

*(**»>)  = 2e  ° -y--  m ,° 


Taking  the  inverse  transform  of  Eg.  (28), 


_ if"  ^ — sin(oD  - cd  ) . . 

- ±-  / ° ° — i — - “1  e*"*  to 

2,t  J.  CO  l®  - ®Q) 


Hence  the  Fourier  components  of  the  incident  pulse  are 


sin(cD  - cd  ) i[cot-HD  — ] 

p(r,t ) = - °1  e ° c 

' * Jt  (cd  - 03  ) 

' 0 ' 


, Sin(o3  - 03  ) - i ( 03-03  )—  ifcDt+03  — ] 

_/ ? 4.  \ 1 O'  v o'c  c 

PI  r,t ) — — ..  — . - - ■ * — 0 0 

' Jt  (03  - 03  ) 

v o' 


In  the  appendix  it  is  proved  that  when  a sinusoid  of  the  form 


P.  = P e 
1 o 


i[“lt+a3l  c] 


is  scattered  from  this  cylinder,  the  scattered  wave  has  the  form 


io3..t  V1  e in  J (o3n  — ) . . 

p = -p  « 1 \ n nv  1 cJ  tt  [2),  r \ 

s p° e L ~rj2 r( — iTT  V (“l  c)  cos  n0 
n=°  V Kc) 


m 


■ • ' 

i 


Thus  if  the  incident  sinusoid  has  the  form  of  Eq.  (29),  the  scattered 
wave  must  have  the  form 


, . ( t i sin(“  - “o'*  ’1(<o-a>oh- 
=)  n (co  - cdJ  6 


oo  . n / a \ *"1 

V enL  JrS^  1T  (2),  r [ ioyt 

/ H v (car)  COS  n9>  e 

U TTsyZaT  n c 


n=o  V (“?) 


(32) 


it 


The  scattered  pulse  is  found  by  integrating  Eq.  (32)  over  all 
frequencies . 


y(?'t)  * s /„ 


2sin(oo  - oo  ) -i(oo-oo  ) — 
' o'  o c 

^ (CO  - O3o)  6 


J '**^2  H (2)(o£)  cos  nej  e^  dno 

H^(«£)  n C'  | 


(33) 


Rewriting  Eq.  (28), 


X(r,a>)  = 2e 


ini  — sin(uwi)  ) 


o c 


(co-a>0) 


By  definition,  the  transfer  function  for  the  scattering  process 
observed  at  r^  is 


H(rt,ci)) 


Y(rt,a^) 

X(rt,a)) 


(35) 


Substituting  the  expressions  for  X(r,cs)  and  Y(r,o))  into  that  for 
H(rt,a>), 


n=o 


(a)—- ) cos  n6. 


(36) 


B.  Example  2:  Transfer  Function  and  Impulse  Response  for  a 

Rigid  Scatterer 

Assume  continuous  s ingle- frequency  scattering  from  a convex 
rigid  body.  The  Kirchoff  Approximation  and  the  physical  optics 
approximation  are  valid.  Let  the  transmitter  and  receiver  be 
coincident  at  the  origin  of  a polar  coordinate  system.  (See  Fig.  3.) 
A point  P has  arbitrary  location  on  the  insonified  surface  2 of  the 
rigid  body,  fi  is  the  unit  normal  to  L at  P,  and  R is  the  distance 
from  0 to  P.  The  smallest  angle  between  OP  and  n is  designated  i(r. 
Assume  that  both  the  transmitter  and  the  receiver  are  omnidirectional 
and  have  a flat  frequency  response.  Designate  as  the  input  of  a 
linear  system  the  electrical  signal  in  the  transmitter.  Let  the 
output  be  the  electrical  signal  in  the  receiver.  Assume  that  the 
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medium  is  one  in  which  the  homogeneous  wave  equation  holds  exactly. 
For  a continuous,  single- frequency  sinusoid  of  unit  amplitude,  the 
input  is  written 


■ mm 


, . \ icot 
s(t)  = e 

This  is  converted  by  the  transducer  into  a spherical  wave  in  the 
water . 


x{r,t)  - i . 

The  scattered  wave  observed  at  the  origin  (i.e.,  the  output  of  the 
system)  is  found  from  the  Kirchoff  integral  solution  to  the  wave 
equation  to  be 


P 

, 2R  1 

r -ln>~r 

e(t)  = - 

1«c  J 

1 

c cos  ifdZ  ) 

£ R^  j 

By  the  process  outlined  in  Example  1,  the  transfer  function  for  this 
system  is 


. r -iciy— 

H(u>)  = Jz  e c cos  i|rd£  . (37) 

With  considerable  manipulation,  the  inverse  Fourier  transform  of  this 
result  yields 

*(t)  = nk  ^5  f H(*  - £2|1  42  > (5®) 

At  JZ  R 
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2R 

where  H(t  - — ) is  the  Heaviside  step  function. 

But  cos  ^ j_s  -the  differential  element  of  a solid  angle  subtending 
R 

the  surface  d£. 

Thus 


h(T)  ' TSc  £2  /„  H<T  ‘ “>®  ’ 

where  M is  the  solid  angle  subtended  at  the  receiver  by  that  portion 
of  the  rigid  body  surface  visible  from  the  receiver. 

The  value  of  the  integral  at  any  time  t,  W(t),  is  the  solid  angle 
subtended  by  that  portion  of  the  visible  surface  between  the  receiver 
and  a spherical  shell  of  radius  (Because  of  the  step  function, 

the  contributions  from  the  other  side  of  the  shell  are  zero.) 

h(f)  = ^“5  W(t)  . (39) 

dx 

If  the  transmitter  and  the  receiver  are  very  far  from 
the  spherical  shell  is  nearly  plane  in  the  region  of 
To  a close  approximation  then, 

w(t)  = ig  A(t)  , 

R 

where  A(t)  is  the  area  of  the  intersection  between  this  plane  and 
the  rigid  body,  and  R is  the  distance  from  the  receiver  to  any  point 
on  the  rigid  body. 

Hence  we  obtain  the  following  result  when  the  transmitter  and  receiver 
are  in  the  far  field. 


(See  Fig.  4.) 
the  scatterer, 
the  scatterer. 
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III.  EXPERIMENTAL  DETERMINATION  OF  THE  TRANSFER 
FUNCTION  AND  IMPULSE  RESPONSE 

As  shewn  in  Chapter  I,  the  impulse  response  and  transfer 
function  of  a linear  system  constitute  a Fourier  transform  pair. 
Consequently,  a determination  of  one  enables  the  other  to  be  calcu- 
lated. The  transfer  function  is,  in  general,  complex  and  may  be 
written  as 

H(cd)  = A(cu)ei^C°^  , 

where 

A(cd)  = |h(cd)  | , and 


q>(cu)  = arg[H(uj)] 


A(cd)  is  the  amplitude  of  the  system  response  to  a unit  sinusoid  at 
frequency  u>  and  cp(ui)  is  the  phase  difference  between  input  and  output. 
A plot  of  H(cd)  vs  cd  can  be  obtained  for  a particular  system  by 
applying  unit  sinusoids  of  various  frequencies  and  measuring  the 
amplitude  and  phase  of  the  outputs . 


If  the  means  were  available,  the  transfer  function  could  be 
obtained  by  applying  concurrently  unit  sinusoids  of  all  frequencies 
(with  linearly  varying  phases),  and  by  measuring  all  output  ampli- 
tudes and  phases . But  because  impulses  cannot  be  physically- 
produced,  this  method  is  impractical.  The  following  method  of 
obtaining  the  transfer  function  involves  considerable  computation. 
A broad  band  pulse  is  scattered  from  the  target.  The  incident  and 
reflected  pulses  are  recorded  on  tape,  digitized,  and  a transfer 
function  is  calculated  by  a computer  from  the  relation 
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Some  pulses  have  spectra  that  are  approximately  flat  over  a 
certain  frequency  range.  (An  example  is  the  large- amplitude  Gaussian 
pulse  with  unit  area,  Eq.  (43).  If  such  a pulse  is  applied  to  a 
linear  system,  the  spectrum  of  the  response  will  be  approximately 
proportional  to  the  transfer  function  in  that  frequency  range.  If 
the  system  responds  mainly  to  inputs  in  the  same  frequency  range, 
a transform  of  the  system  function  thus  derived  will  approximate 
the  impulse  response. 

Another  method  of  finding  the  impulse  response  is  based  on 
Eq.  (19).  For  an  input  x(t)  (~  white  noise)  and  an  output  y(t), 

h(t)  = aRyx(t)  , 

where 

a = constant 

If  the  input  approximates  white  noise,  the  output  will  have  a 
cross  correlation  roughly  proportional  to  the  impulse  response. 


IV . TARGET  STRENGTH  AND  THE  TRANSFER  FUNCTION 


One  goal  of  scattering  studies  is  to  associate  with  each  target 
a set  of  numbers  characteristic  of  that  target's  scattering  behavior. 
One  such  set  of  numbers  is  target  strength  measured  for  various 
aspects  and  frequencies.  Target  strength  is  defined  as  the  ratio 
of  the  reflected  intensity  referred  to  1 yd  from  the  effective 
acoustic  center  to  incident  acoustic  intensity  (assuming  plane  waves). 
It  can  be  measured  as  follows . 

Let  a source  at  Q,  a distance  d from  a target  S,  send  out 
cw  plane  waves  perceived  at  a hydrophone  R.  The  hydrophone  lies  at 
a distance  r from  the  target.  (See  Fig.  6.)  Let  Pq  be  the  ampli- 
tude of  the  steady- state  part  of  the  incident  pulse;  and  let  be 

the  amplitude  of  the  steady- state  part  of  the  scattered  pulse 
-¥  “► 
measured  at  r.  The  target  strength  for  a particular  r and  o>  will  be 

P- 

T^  = 20  logl0  + 20  log(r  - l)  . (4l) 

o 

The  transfer  function  can  be  obtained  by  methods  discussed  in 
Chapter  III.  No  matter  how  it  is  determined,  the  transfer  function 
can  be  written  in  the  form 


«(?,»)  = ps  , (w) 

o 

where  P*  and  P are  the  quantities  described  above, 
o 

From  Eqs . (4l)  and  (42)  it  is  apparent  that  the  target  strength  is 
the  amplitude  part  of  the  transfer  function  corrected  for  distance 
and  expressed  in  decibels.  It  does  not  contain  phase  information. 
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V.  EXPLOSIVE  PULSES 


The  unit  impulse  or  Dirac  delta  "function"  does  not  represent  a 
physically  realizable  waveform.  Since  the  impulse  response  of  a 
linear  system  is  the  output  of  that  system  corresponding  to  a unit 
impulse  input,  a physical  system  can  never  have  its  own  impulse 
response  as  an  output.  Because  all  physical  systems  behave  non- 
linearly  when  large  amplitude  inputs  are  employed,  the  impulse 
response  can  not  even  be  considered  as  a limiting  case  of  the  output 
of  a physical  system. 


The  impulse  response  is  a coded  statement  about  the  response 
of  a physical  system  to  inputs  of  arbitrary  waveform.  It  may  be 
decoded  and  applied  to  specific  cases  using  linear  system  theory. 
Since  the  impulse  response  is  characteristic  of  a target's  acoustic 
behavior,  it  might  be  assumed  that  the  response  to  explosive  pulses 
would  also  be  characteristic.  Strangely  enough  this  appears  to  be 
the  case.  It  will  be  shown  that  for  at  least  one  short-duration, 
large- amplitude  input,  the  output  of  a soft  cylinder  resembles  its 
impulse  response. 


Suppose  the  Gaussian  pulse  of  unit  area 


-i/2(S& *)2 

e / a 


(*5) 


is  incident  on  the  cylinder  in  Fig.  2.  The  scattered  pulse  must  be 


y(r,t)  = F_1  [X(r,u>)  H(r,u>)] 
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The  Fourier  transform  of  the  incident  pulse  was  found  to  be 
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Now  because 
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the  output  must  be 
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But  the  impulse  response  is  the  inverse  transform  of  H(r,a>), 
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If*  the  input  is  of  large  emplitude,  o must  be  a small  number. 

This  means  that  for  some  frequency  range,  e _1/2^ CT/C  will  have  a 

value  close  to  1.  If  the  target  or  receiver  is  band-limited  to  the 

same  range  Acs,  e-1/2^a/c  ^ can  be  replaced  by  1 in  Eq.  (47).  The 
resulting  expression  differs  from  the  impulse  response  only  by  an 

, a1 

- ixj>~ 

amplitude  scale  factor  l/c  and  a time  delay  factor  e c (depending 
on  the  location  of  the  observer).  Implicit  in  this  result  is  the 
assumption  that  the  large  amplitude  of  the  input  does  not  result  in 
significant  nonlinear  effects  on  the  scattering  problem.  Based  on 
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these  assumptions,  the  scattered  wave  will  resemble  the  impulse 
response  of  the  target  (as  seen  by  the  receiver)  for  at  least  one 
input  waveshape  of  large  amplitude  and  short  duration. 

This  similarity  and  other  factors  have  led  to  experimental 
investigations  of  the  propagation  and  scattering  of  explosive  pulses. 
Work  with  one  impulsive  source  (see  Cohen  and  Winder,  1966)  indi- 
cates high  reproducibility  of  waveform  at  short  range  (less  than 
100  yd).  For  the  same  source,  TO  percent  of  the  energy  is  contained 
in  the  1 to  20  kHz  band.  As  range  increases,  the  high-frequency 
content  decreases.  Distortion  of  the  pulse  in  transmission  is  con- 
siderable. Fortunately  the  distortion  appears  reproducible  for 
short  ranges,  and  the  pulse  retains  its  impulsive  character  out  to 
at  least  100  yd.  The  beamwidth  of  the  peak  pressure  pattern  at  the 
-3  dB  level  was  about  30  deg  at  a range  of  1 yd. 


Consider  a soft  cylinder  of  infinite  length  and  radius  a, 
symmetric  about  the  7-axis.  (See  Fig.  J.)  A plane  sinusoid  of 
frequency  co  is  incident  from  the  right. 


i(cot+kr  cos 


The  net  pressure  in  the  medium  is 


where  P is  the  contribution  from  the  scattered  wave,  and  P must 

b S 

satisfy  the  wave  equation;  i.e., 


The  solution  is 


The  boundary  conditions  to  be  satisfied  are 
l)  The  cylinder  is  soft. 


2)  lim  P is  a plane  type  wave  traveling  directly  away  from 
s 

X*-K>0 

the  7- axis. 


Considering  the  second  boundary  condition  first,  in  view  of  Eq.  (Ak), 
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Using  asymptotic  expansions  for  and  Y^,  and  assuming  r > n,  Eq.  (A6) 
can  be  rewritten 
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Eq.  (A7) 


+ rin  sin(kr  - rr  ~ tj)]  cos  n©J 

If  this  expression  is  to  represent  a plane  type  wave,  we  must  have 
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For  an  outward  traveling  wave. 


tj  = -i| 
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Accordingly,  Eq.  (A1)-)  can  be  rewritten 
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Eq.  (A8) 
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The  incident  wave  can  be  expanded  into  the  following  series  of  Bessel 
functions 


P.  = P e*^  / € in  J (kr)  cos  n 6 , 
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where 

en  = 1»  n = 0 > 

e = 2,  n > 0 
n 7 

Consider  the  first  boundary  condition 


P.  + P = 0@r  = a 
i s 


From  Eqs . (A9),  (A10),  and  (A5),  it  must  be  true  that 
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Eq.  (A10) 


Eq.  (A5) 


Eq.  (All) 


This  will  be  true,  in  general,  only  if 


V.1*  Jn(ka) 

Hnl2,(ka) 


Eq.  (A12) 
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Substituting  Eq.  (A12)  into  Eq.  (A9)  results  in  the  solution 


(The  reverse  of  this  page  is  blank.) 
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